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Abstract 

Although optical metamaterials that show artificial magnetism are mesoscopic systems, they are 
frequently described in terms of effective material parameters. But due to intrinsic nonlocal (or 
spatially dispersive) effects it may be anticipated that this approach is usually only a crude approx- 
imation and is physically meaningless. In order to study the limitations regarding the assignment 
of effective material parameters, we introduce a technique to retrieve the frequency-dependent el- 
ements of the effective permittivity and permeability tensors for arbitrary angles of incidence and 
apply the method exemplarily to the fishnet metamaterial. It turns out that for the fishnet meta- 
material, genuine effective material parameters can only be introduced if quite stringent constraints 
are imposed on the wavelength/unit cell size ratio. Unfortunately they are only met far away from 
the resonances that induce a magnetic response required for many envisioned applications of such 
a fishnet metamaterial. Our work clearly indicates that the mesoscopic nature and the related 
spatial dispersion of contemporary optical metamaterials that show artificial magnetism prohibits 
the meaningful introduction of conventional effective material parameters. 

PACS numbers: 78.20.Bh, 78.20.Ci, 41.20.Jb 
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The introduction of nanostructured metamaterials (MMs) into optics potentially opens 
the door to a fairly comprehensive control of light propagation. During the past several years 
much effort has been devoted to achieve this goal and two major research fields may be dis- 
tinguished. On the one hand, advances in nanotechnology provide ever smaller and more 
complex structures which constitute quite involved nanostructured media. On the other 
hand, optics in media with unprecedented effective material equations has been investigated 
purely theoretically and novel effects have been revealed^. The desirable assignment of ef- 
fective material parameters to a specific MM would bridge the gap between both approaches 
and allow to link a fabricated structure to a particular effective constitutive relation^-. Fur- 
thermore it would appreciably facilitate the description of light propagation in optical MMs 
and their combination with other optical materials because canonical Maxwell boundary 
conditions could be applied. If we wish to understand under a MM an artificial medium 
made of periodically or non-periodically arranged metaatoms which allows to control the 
properties of light propagation predominantly by the chosen geometry of the metaatoms, 
many different MMs can be envisioned which all have peculiar aspects if effective properties 
shall be designed^ 7 . To avoid any misunderstanding, we wish to restrict our considerations 
in the following on metamaterials that show the effect of an artificial magnetism and which 
shall operate at optical frequencies. Such property is often at the focus of interest since the 
media would enable optical phenomena that contradict our common perception of how light 
propagates. 

For such special type of MM, there is, however, a serious issue which might prevent this 
simplified description in terms of effective material parameters. In general, typical optical 
MMs are mesoscopic where the vacuum wavelength is only a few times larger than the unit 
cell. In such systems the optical response may be reasonably described by induced currents, 
which nonlocally depend on the electric field. In Fourier space this leads to a spatially 
dispersive conductivity, as discussed in the work of Serdyukov et al. [8]. At this stage 
it is not required to distinguish between polarization (~ JjP) an d magnetization currents 
(~ V x M), however, this becomes important if either of them becomes resonant in the 
nanostructure. Provided that this spatial dispersion is weak, the constitutive relation be- 
tween j and E can be expanded up to the second order. Since the fields D and H cannot 
be defined uniquely, spatial derivatives of E may be replaced in favor of B. As a result, 
two constitutive relations D(E, B) and H(E, B) emerge with tensorial, but only frequency 
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dependent coefficients^. First order terms lead to magnetoelectric (E <=± H) coupling (bian- 
isotropy, chirality). Second order terms to anisotropic, but spatially nondispersive relations 
between both D and E and H and B. To sum up, these so-called bianisotropic constitutive 
relations are the most general ones for a weak spatially dispersive conductivity in a nanos- 
tructured material. From a physical point of view it is appealing that the nonlocal relation 
between the electric field and the induced currents is the very source for the effective chiral 
and magnetic (fl(u)) properties of MMs. From a technical point of view, with these spatially 
nondispersive constitutive relations at hand, standard boundary conditions^ can be used to 
solve macroscopic Maxwell's equations in layered media. This has a big advantage compared 
to the rather involved procedure for spatially dispersive constitutive relations which require 
the use of so-called additionally boundary conditions 10 . In mirror-symmetric (nonchiral) 
media first order terms in the expansion vanish and the magneto-electric coupling disap- 
pears. The MM may then be described by two material tensors i(u) and fi(oo). 
Here we aim at introducing a simple criterion that tells us if this condition is fulfilled. 
In general, one has to develop an approach to retrieve these tensor elements from reflec- 
tion/transmission data. However, it will turn out that the very calculation of these param- 
eters is not required. Only if the criterion is fulfilled one has to proceed with the retrieval 
algorithm to calculate effective parameters which are then independent of the incidence angle 
and may be termed effective material parameters. On the other hand, if the criterion is not 
fulfilled the assignment of an effective permittivity s(oj) and permeability fi(u) is pointless. 
This means physically that the assumption of weak spatial dispersion is violated and the 
effective material parameters would become spatially dispersive making the approach used 
inconsistent. Thus the aim of this work is not the retrieval of parameters but to evaluate if 
a certain MM may be described by effective material parameters. 

But in any case the recently introduced retrieval approach for isotropic materials at arbi- 
trary incidence^ has to be generalized towards anisotropic media. The main advantage of 
the present approach is that all tensor elements can be determined without requiring ex- 
plicitly that the propagation direction coincides with a crystallographic axis. Hence, the 
approach may be applied to all currently fabricated MMs. 

To start with, we assume that the metallic structures (split ring, fishnet, etc.) involved are 
reciprocal and not intrinsically magnetic (/^ me tai = !)• Their response to the electromagnetic 
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field can be completely described by a current nonlocally induced by the electric field as^ 

j(r,w) = J R(r,r',uj)E(r',cj)dr', (1) 



v 



where the dyadic R(r, r', u) describes the nonlocal response of the medium. However this 
approach is not practical. For weak nonlocality (or spatial dispersion) one may rather expand 
Eq.flTJ up to the second order. For improving the readability of the paper we provide our 
derivation in following the lines in Ref. bi We obtain 



dEi dE l 

O-klEi + Ofci m — h Ckh 



(2) 
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where Einstein notation is applied and bki m = — hkm- The factor icj has been introduced 
for convenience, so that the terms in the brackets represent the induced polarization. Now 
the constitutive relations read as 

, \ / r \ en dE\ <9 2 Ei , . 

J^k\ v i = (e 0kl + Clkl) Ei + Oklm-R 1" CklmnJ, > V>) 

ox m ox m ox n 
H= B//i . (4) 

The first term leads to the anisotropic permittivity Ski = Ski + a>ki/ £ o- The second term 
accounts for magnetoelectric coupling and vanishes for media with three orthogonal planes 
of mirror symmetry, i.e. media that are purely anisotropic in the quasi-static limit. To 
proceed we require the coefficients Cki mn to obey the following relation: 

d 2 E i 

Cklmn -KZ ' = [V X (7V X E)] fc (5) 
UX rn UX n 

where the components 7^ are related to the coefficients Cki mn - Since Cki mn = Cikmn holds, 
the tensor 7 is symmetric (7^ = 7^). 
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Because Maxwell's equations are invariant with respect to the transformations 
D + V x Q, H' = H — iuQ we can rewrite Eqs. [3] and H] in using Q(r, lo) = —7V x E = iuj~B 
to obtain the ultimate constitutive relations 

D(r, u) = e i(uj)E(r, u), B(r, u) = yU /i(w)H(r, u). (6) 

These equations represent our point of departure. They reflect that a mesoscopic metallic 
structure with a weak nonlocal response (weak spatial dispersion) can be likewise treated as 
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an effective homogeneous, anisotropic but magnetic medium where the magnetic properties 
merely originate from the nonlocal response R(r,r',u): 

m = (Ski ~ /ioW 2 7fci) _1 - ( 7 ) 

Only if the constitutive relations of a medium without bianisotropy can be cast in the form 
of Eqs. (Q, the usual boundary conditions are applicable where the tangential components 
E t and H t as well as the normal components D n and B n are continuous. 
The strategy of our work is as follows: We develop a retrieval algorithm for anisotropic, ho- 
mogeneous and local media which relies on the rigorously calculated reflection/transmission 
data from a MM slab for transverse wave vectors that extend even into the evanescent 
domain. Within this algorithm we identify a quantity a(u) which contains the essential 
information. If this quantity depends only on frequency and not on the wave vector, the 
effective permittivity and permeability tensor elements, not calculated at this step, will 
exhibit the same feature and can be calculated in a further step and assigned to the meta- 
material. If this is not the case the ratio unit cell size/wavelength has to be decreased until 
this criterion is met. Here we restrict ourselves to media where both material tensors can 
be simultaneously diagonalized as 

e = di&g{e x , e y , e z }, fx = diag{/i x , /j y , fx z }. (8) 

We align the coordinate system and therefore all interfaces to the crystallographic axis 




FIG. 1: (color online) a) Schematic view of the single fishnet layer together with the four principal 
directions for the retrieval, b) Unit cell of the fishnet with P x = P y = 600 nm, W x = 284 nm, 
W y = 500 nm embedded in air. The thicknesses of the silver and the intermediate MgF 2 (n = 1.38) 
layer are d,Ag = 45 nm and d = 30 nm, respectively. 

of the effective anisotropic medium. The incident light shall consist of monochromatic 
plane waves whose wavevector is perpendicular to at least one coordinate axis. Then the 
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eigenmodes of the medium can be decomposed into decoupled TE and TM modes and the 
reflection/transmission problem is equivalent to that of an isotropic medium. The only 
difference is the different propagation constant for each eigenmode. For a particular eigen- 
polarization the reflection and transmission coefficients in terms of the normal wavevector 
component have the same form as in the isotropic case. For varying the incidence plane and 
the polarization only certain quantities have to be exchanged as indicated in Tab. [H The 
transmission and reflection coefficients for the electric field read as 

^ = t(k s + k c ) cos(kd) - + k s k c ) sin(kd) ^ 

„/, c , _ £(k s - k c ) cos(kd) + i(^ 2 - k s k c ) sin(fcd) 
1 ' U ~ £(k s + k c ) cos(kd) - i(e + k s k c ) sin(kd) [ J 

where the following abbreviations have been used: 



k s , c = a s > c k s z ' c , k = H, A TE = 1, A TM = J £ ^ (11) 



with £ = a f k f z , where k\ = J - k 2 t ■ 7* (12) 

is the normal component of the wavevector in medium '%' and k t = {k x k y ) is its conserved 
tangential component. The superscripts i G {s, f, c} denote substrate, film and cladding. In 
what follows we assume that substrate and cladding are isotropic and nonmagnetic where 
a TE = 1; «tm = l/£ s ' c > P s ' c = £S,C an d Y' c = 1- F° r the sake of clarity we drop the superscript 
'/' and write a^(u) = a(u>), 0^(oo) = f3(u) and 7^(w) = 7(0;). These coefficients are 
related to different combinations of tensor components of the permittivity and permeability 
depending on the polarization and the incidence plane, see Table HI For the sake of brevity 
the frequency dependence will be kept in mind but not explicitly written in the following. 

The effective material slab is then fully characterized by the parameters k and £. Note 
that throughout the paper k is the normal component of the wavevector in the slab. By 
inverting Eqs. ([9]) and ( flQ|) one obtains 

kd = iarccos 7— . „. ri . — —. — T + Lvmx (Id) 

\{T/A)[k s {l-R) + k c {l + R))J 1 ; 

with m G Z and 



j k 2 s {R-iy-kl{T/AY 
K ~ V {R+lf-{T/Af ' K } 
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TABLE I: Substitution table for the relevant coefficients depending on the polarization and the 
incidence plane. 

The sign of k and £ and the branch order m are determined by the usual physical 
constraints^. The quantities k and £ can be uniquely determined and are the final ef- 
fective wave parameters. They are independent of the thickness d of the slab, provided that 
they already converged towards the bulk datar^. Then these wave parameters must coin- 
cide with those provided by the dispersion relation of the fundamental Bloch mode^. It is 
evident that these effective wave parameters describe properties of the fundamental Bloch 
mode in the infinite lattice formed by periodically arranging a single fishnet layer— They 
still depend on the propagation direction, the angle of incidence and the polarization state 
as in any anisotropic medium. 

Now the criterion central to this work can be formulated. Since a (a;) = £/k (see Eq. [T2l) 
is related to the effective material parameter tensor components (see Tab. UJ) it has to be 
independent of the angle of incidence although £ and k will strongly depend on it. Actually 
this is a very simple criterion that we can use to evaluate the validity of the effective ma- 
terial approach just by calculating £ and k from the transmission/reflection data by using 
Eqs. (d3D and (ED. 

Only in the case that the above criterion is fulfilled the very parameter retrieval can be 
performed, otherwise their assignment is pointless. This parameter retrieval may then be 
performed as follows. If the wave parameters k and £ are determined for a certain polariza- 
tion and incidence plane, the parameter a is given by a = £/k. The parameter /3 follows 
then from Eq. ( JT2l) at normal incidence which reads as k = —yf$. With /3 being constant, 
the remaining parameter 7 can be retrieved as a function of kt from Eq. ( IT2l) as 7 = — p — . 
Now having a, (3, and 7, the classification in Tab. Ucan be used to retrieve the effective ma- 
terial parameters for the respective polarization and incidence plane. It can be recognized 
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that each parameter can be retrieved twice; providing a possibility to double-check results. 
If both results agree the model of weak spatial dispersion holds and the effective material 
parameters are meaningful. 

The outlined procedure is exemplarily applied to the fishnet MM. Geometrical parameters 
are given in Fig. [TJ The metal layers are assumed to be made of silver—. To identify the 
limits of the effective parameter description, we varied the wavelength to unit cell size ratio 
by scaling the structure and calculated a = £/k. Hence, all parameters defined in Fig. [1] 
are scaled by a factor /. Results of the numerical simulations are summarized in Fig. [2] 
where the real part of the propagation constant k is shown as a function of the wavelength 
and the scaling factor for normal incidence. The dark-colored area exhibits a negative prop- 
agation constant and can thus be considered magnetically active (antisymmetric plasmon 
resonance). In the light-colored area the fishnet exhibits a plasma-like effective permittivity 
with negligible dispersion in the permeability. The green solid line indicates the region of 
5f(/x) < at normal incidence, i.e. the region of double negativity. It is evident that the 

a) b) 



FIG. 2: (color online) a) Real part of the propagation constant k (in /xm _1 ) of the fishnet structure 
for normal incidence. Additionally the lines of constant relative deviation of a are shown (black 
cross - 0.1%, cyan asterisk - 5%). The relative deviation is calculated for varying kt = k y in TE 
polarization (a = \i~ ) being the preferential operating polarization state for the fishnet. The 
green solid line indicates the region of K(^) < 0. b) zoomed domain of interest from a). 

smaller the structure the smaller the resonance wavelength and the weaker the resonance 
strength. For very small scaling factors (/ < 0.1, i.e. unit cell size of less than 60 nm) where 
the quasi-static limit is reached, the resonance is almost wavelength independent but also 
tends to disappear; being in agreement with findings for split-rings— 1 ^. 
Having identified the area where resonances are occurring it is now interesting to disclose 
where the effective parameter description (a must be invariant) may be applied. To this 
end we calculated for discrete wavelengths the parameter a depending on the transverse 
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wavevector (k t = 0..1.2/co) including at least partially the evanescent spectrum. The relative 
deviation (max \a(kt) — a(0)|)/|a(0)|) serves as a measure to characterize the variation of 
a and is displayed in FigfSJ Two bounds for this deviation are considered, 0.1% (almost 
ideal assignment of effective material parameters possible) and 5% (assignment of effective 
parameters might be still feasible). Clearly, close to the resonance this deviation is strongest 
for a fixed scaling factor. It is evident that for a required deviation of 0.1% effective material 
parameters can be only introduced when there are no magnetic resonances (no lefthandness). 
Since an effective description holds only for almost invariant a this deviation should be as 
small as possible but values of 5% might be tolerable at most. This condition requires a 
scaling factor of about 0.15 in the resonance region resulting in a wavelength to cell size 
ratio of \/P > 10. Evidently, the resonance strength is very weak in this domain leading 
to a non-magnetic response of the material (lefthandness occurs only because of the large 
imaginary parts). This is consistent with the assumption that such small structures can 
be described in the quasi-static limit where no magnetic response is observed. Hence, the 
result of our studies is quite discouraging, namely: a sufficiently strong magnetic response 
(5f(/x) < 0) requires a certain minimum unit cell size/wavelength ratio (about 1:4 in case of 
the fishnet), but this mesoscopic structure must not be described by conventional frequency- 
dependent effective permittivity and permeability tensors. We have proven this for a fishnet 
structure, but since all present optical MMs rely on similar resonances we conclude that this 
tendency may hold in general. 

To sum up, based on the assumption of a weakly spatially dispersive conductivity in MM 
unit cells we have developed a method to verify /falsify an effective anisotropic medium 
description of MMs. We have shown that a prototypical magnetically active, and thus po- 
tentially negative index, material, namely the fishnet, cannot be described as a homogeneous 
anisotropic medium in the relevant resonance region. By varying the ratio of wavelength 
to cell size we have elaborated the limitations of the weak spatial dispersion assumption. 
There is a trade off: If the spatial dispersion is weak and the material parameters have the 
usual meaning, the antisymmetric plasmonic resonance, which is responsible for magnetic 
activity, is also weak or disappears. Our work clearly indicates that for optical MMs the 
commonly assigned effective parameters do not have the physical meaning of conventional 
material parameters. 

This work was partially supported by the German Federal Ministry of Education and Re- 
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